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Boundary layer stability was analyzed for the HIFiRE-1 flight vehicle geometry for ground tests con- 
ducted at the CUBRC LENS I hypersonic shock test facility and the Langley Research Center (LaRC) 20- 
inch Mach 6 Tunnel. Boundary layer stability results were compared to transition onset location obtained 
from discrete heat transfer measurements from thin film gauges during the CUBRC test and spatially con- 
tinuous heat transfer measurements from thermal phosphor paint data during the LaRC test. The focus of 
this analysis was on conditions at non-zero angles of attack as stability analysis has already been performed 
at zero degrees angle of attack. Also, the transition onset data obtained during flight testing was at non- 
zero angles of attack, so this analysis could be expanded in the future to include the results of the flight test 
data. Stability analysis was performed using the 2D parabolized stability software suite STABL (Stability 
and Transition Analysis for Hypersonic Boundary Layers) developed at the University of Minnesota and 
the mean flow solutions were computed using the DPLR flnite volume Navier-Stokes computational fluid 
dynamics (CFD) solver. A center line slice of the 3D mean flow solution was used for the stability analysis 
to incorporate the angle of attack effects while still taking advantage of the 2D STABL software suite. The 
N-factors at transition onset and the value of Reg /Me, commonly used to predict boundary layer transition 
onset, were compared for all conditions analyzed. 

Ground test data was analyzed at Mach 7.2 and Mach 6.0 and angles of attack of 1°, 3° and 5°. At these 
conditions, the flow was found to be second mode dominant for the HIFiRE-1 slender cone geometry. On the 
leeward side of the vehicle, a strong trend of transition onset location with angle of attack was observed as 
the boundary layer on the leeward side of the vehicle developed inflection points at streamwise positions on 
the vehicle that correlated to angle of attack. Inflection points are a strong instability mechanism that lead 
to rapid breakdown and transition to turbulence. The transition onset location on the windward side of the 
vehicle displayed no trend with angle of attack or freestream Reynolds number and transition was observed 
farther down the vehicle than observed on the leeward side of the vehicle. In analysis of both windward 
and leeward sides of the vehicle, use of the N factor methodology to develop trends to predict boundary 
layer transition onset showed improvements over the Reg /M^ empirical correlation methodology. Stronger 
correlations and less scatter in the data were observed when using the N factor method for these cases. 

I. Introduction 

The HIFiRE-1 test program was developed to advanee researeh in weapons systems by demonstrating instrumen- 
tation capabilities on a flight vehicle ^ . Laminar to turbulent transition and turbulent separated shock-boundary layer 
interaction were the flight objectives chosen to demonstrate those capabilities and the HIFiRE-1 flight geometry was 
designed with these test objectives in mind. First, a 7° half-angle cone with a 2.5 mm nose radius was selected for 
the fore body and instrumented to observe natural and tripped transition to turbulence. At the aft end of the cone 
was a straight, cylindrical section to allow for the turbulent flow to separate, and then a 33° flare was attached at the 
end of the cylinder and was densely instrumented to obtain data on the resulting turbulent separated shock-boundary 
layer interaction region. The flight geometry can be seen in Figure 1. 
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A. Review of Ground Testing 

The HIFiRE-1 geometry was tested in two ground test faeilities. The CUBRC LENS I hypersonie shoek test faeility 
in Buffalo, New York tested the full seale geometry at flight Reynolds numbers to duplieate the reentry trajeetory.^ 
The Langley Researeh Center tested a sealed flight geometry in the 20-ineh Maeh 6 Tunnel. ^ Information on both 
of these tests ean be found in other publieations in more detail, so only a brief overview will be eovered here. 

During the test at the CUBRC LENS I faeility, two entries were eompleted with varying nose radii and flare 
angles to optimize the geometry to meet the flight test objeetives. Angles of attaek of 0°, 1°, and 5° were also 
tested. Diserete thin film sensors were plaeed at four azimuthal loeations to eapture the transition behavior at angles 
of attaek. Runs seleeted from the CUBRC test series for this analysis eonsisted only of runs where the geometry 
eonfiguration matehed the eonfiguration flown during the flight test. 

The LaRC test series in the 20-ineh Maeh 6 Tunnel tested two models eaeh with two different nose radii (four 
models in total). The first was a 20% seale model of the entire flight vehiele and the seeond was a 35% seale 7° 
half-angle eone without the eylinder and flare. Thermal phosphor paint was used to eapture spatially eontinuous 
heat transfer data to gain global transition behavior at 0°, ±3°, and ±5° angles of attaek. The runs at positive and 
negative non-zero angles of attaek at eaeh run eondition were designed to eapture transition data on both windward 
and leeward sides of the vehiele. 

While many runs were eompleted for the tests eondueted at CUBRC and LaRC to meet different test objeetives 
and to gather transition data for a variety of eonditions. Table 1 lists the runs and eonditions ehosen for the analysis 
in this publieation. These eonditions were ehosen beeause data was gathered at non-zero angles of attaek where 
transition was observed. Runs were ehosen to get transition information on both the windward and leeward sides of 
the vehiele, however, not all eases transitioned on both sides of the vehiele. 


Table 1. Ground Test Data Analyzed 


Lacility 

Run 

Mach 

Re/m [xlO-6] 

a[° 

CUBRC 

10 

7.152 

10.23 

1 

CUBRC 

17 

7.140 

10.43 

1 

LaRC 

10/22 

5.977 

8.11 

3 

LaRC 

12/24 

6.007 

13.38 

3 

LaRC 

63/65 

6.028 

18.68 

3 

LaRC 

16/28 

5.977 

8.02 

5 

LaRC 

19/29 

6.006 

13.03 

5 

LaRC 

64/30 

6.028 

18.63 

5 

CUBRC 

11 

7.176 

10.65 

5 

CUBRC 

13 

7.155 

10.37 

5 


Also of importanee is the nose radius of eaeh of the models beeause the nose radius plays a signifieant role in 
the stability of the boundary layer. The model used for the CUBRC test was at the flight seale, so the nose radius 
of 2.5 mm was the same for the test artiele in the tunnel as well as the flight vehiele. While the LaRC test series 
used a sealed model, the nose radius was not direetly sealed the same amount. The LaRC tests used models with a 
1.19 mm nose radius and models with a 2.1 mm nose radius for their testing. When sealed up to the flight geometry, 
this would result in a blunter nose tip than the CUBRC test artiele and the flight geometry. A blunt nose tip serves 
to stabilize the boundary layer and delay transition in most eases and must be taken into aeeount in the analysis 
presented in this paper. 

B. Empirical Transition Prediction Methods: Correlation and Physics-Based Analyses 

Predicting boundary layer transition is a difficult problem that many researchers have spent careers trying to solve. 
Boundary layer transition is a complex process that can be governed by fluctuations in freestream properties, surface 
roughness, surface curvature, boundary layer properties and many other geometric and flow phenomena. While 
many advances have been made in determining the characteristics of a boundary layer and its stability, none have the 
capability to predict transition as a stand-alone method. ^ Even more advanced physics-based methods still require 

2 of 12 


American Institute of Aeronautics and Astronautics 



transition onset information from testing. Depending on the applieation, a strong ease ean be made for empirieal 
eorrelations based on boundary layer properties and equally strong eases ean be made for physies-based methods of 
analyzing boundary layer stability to prediet transition. Empirieal eorrelations based on boundary layer properties 
are very useful in engineering level vehiele design or performanee evaluations when simple and effeetive tools 
are needed to assess a design or situation quiekly^, while physies-based methods are useful for understanding the 
meehanisms that eause boundary layer transition as well as how the different modes internet with eaeh other or the 
amplitude of disturbanees grows or deeays. 

One popular method for predieting the onset of boundary layer transition empirieally eorrelates transition onset 
loeation observed in test data to boundary layer properties at that loeation from a eomputational simulation. However, 
empirieal eorrelations based on boundary layer properties are only valid for similar geometries at similar eonditions 
and usually eannot be generalized to multiple different geometries or test eonditions.^ Two popular boundary layer 
properties used to prediet the onset of boundary layer transition on smooth bodies are Reg, the Reynolds number 
based on the momentum thiekness of the boundary layer, and Reg /Me, the momentum thiekness Reynolds number 
divided by the Maeh number at the edge of the boundary layer. Transition onset loeation ean also be eorrelated to 
many other boundary layer properties, some of whieh are addressed by Campbell, et. al.^, but only Reg /Me will be 
addressed in the work presented here. 

While eorrelations based on boundary layer properties ean be useful for engineering problems, they do not eap- 
ture mueh of the physies of the boundary layer and they do not take into eonsideration many of the variables affeeting 
the growth of disturbanees in the boundary layer whieh ean lead to transition onset. Physies-based methods sueh as 
the method and the parabolized stability equations (PSE) have the ability to ineorporate the physies governing 
the mean flow to traek the growth of disturbanees within the boundary layer that may lead to transition. A higher- 
fidelity meehanism based approaeh to boundary layer stability eould provide more aeeurate methods for predieting 
boundary layer transition onset. A better understanding of the transition meehanisms, however, is required to vali- 
date these methods. Many eonventional wind tunnels produee transition data that is highly influeneed by the noise 
propagated from the turbulent boundary layers that develop on the tunnel walls. ^ As a result, transition data obtained 
in quiet tunnels where the boundary layer on the nozzle walls remains laminar and radiated noise and freestream 
fluetuations are greatly redueed provides validation data for transition from other instability meehanisms. Measuring 
those instability meehanisms presents its own set of ehallenges, though. Reshotko ^ outlines a set of guidelines than 
ean be used for setting up experiments to gather transition data. More and better data on the meehanisms that eause 
transition from a wide variety of sourees would be very benefieial to advaneing boundary layer stability predietions 
and eould be used to validate and inform the physies-based models eurrently used. 

It should be emphasized at this point that it is not the purpose of this publieation to say whieh method is better 
for all applieations. This paper addresses two methods on a single geometry and it is not assumed that the results 
ean aeeurately be extrapolated to another geometry or other eonditions. 

II. Methodology 

All ground test eonditions seleeted for this analysis were at a non-zero angles of attaek to examine the effeets of 
angle of attaek on the boundary layer stability and transition onset loeation. The boundary layer transition data ob- 
tained during the HIEiRE-1 flight test was also at non- zero angle of attaek for both aseent and deseent measurements, 
so future analysis of the flight data eould inelude this methodology. 

The stability analysis was eompleted using the stability modeling software, STABE (Stability and Transition 
Analysis for Hypersonie Boundary Eayers), developed at the University of Minnesota, whieh uses the linear PSE 
(parabolized stability equation) method in PSE-Chem. ^ This method uses a laminar basie state solution and assumes 
the small perturbations in flow properties of the form of Equation 1 for pressure, density, temperature, and three 
eomponents of veloeity. 
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( 1 ) 


p = p + p 
p = p + p' 

T = f + T' 

u = U + u' 

The perturbed variables are then substituted into the boundary layer equations and the system of equations is 
parabolized to obtain the Parabolized Stability Equations (PSE), which reduce to Equation 2 for the two-dimensional 
case of importance here. 


(i- + eL'^ q + eMi — — r (2) 

In the above equation, q is the vector of flow variables, L, L', and M\ are operators in the off-body direction, ^ 
allows for slow variations in the streamwise direction, and r on the right hand side includes the non-linear terms. 
The eigenvalues of the system are found to determine the real and complex disturbance frequencies that are used 
to calculate the amplification of each disturbance. A detailed derivation of the PSE methodology can be found in 
Parabolized Stability Equations by Thorwald Herbert^. 

The growth and decay of those imposed perturbations are then tracked down the length of the body as a function 
of disturbance frequency. The amplitude of each disturbance is integrated using Equation 3 to determine the N factor 
at each point along the body and a maximum N factor curve is created. Eigure 2 shows a sample maximum N factor 
curve. 


r IdA ^ ,A 
= I dx = In—-- 

A dx Aq 


'XO 


(3) 


N Factor Curve for Each Disturbance Frequency 
Maximum N Factor Curve 



x[m] 


Figure 2. Maximum N Factor Curve 

Mean flow solutions were computed using the DPER finite volume Navier-Stokes flow solver on a grid with 
?s8.5 million grid cells. Viscosity was modeled using Sutherland’s Eaw and thermal conductivity was computed 
assuming a constant Prandlt Number. Eor simplicity, the aft end of the geometry was excluded from the analysis and 
only a cone was simulated. The three-dimensional nature of the flow field for a cone at angle of attack can be seen in 
the surface heat transfer and surface shear lines in Eigure 3. On the leeward side of the vehicle, two vortices develop 
as the flow diverges from the windward side. The vortices pull high momentum fluid from high in the boundary 
layer down towards the surface of the vehicle which modifies the mean flow solution and introduces instabilities into 
the boundary layer. This must be noted as it will be important in later analysis. 
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(a) Surface heat transfer 


(b) Surface shear lines 


Figure 3. Solution at 5° angle of attack. 


The STABL software used for analyzing the boundary layer stability is a 2D code, but a cone at a non- zero angle 
of attack produces a 3D flow field. For this reason, a centerline slice was extracted from a 3D mean flow solution 
for the stability analysis of each condition. The centerline slice contained 157 cells in the off-body direction and 
approximately 75-100 of those cells resolved the boundary layer. The grid was split at the stagnation point to allow 
for analysis of the windward and leeward stability separately. The flow on the centerline pitch plane used for these 
analyses can be seen in Figure 4. The frequencies of the disturbances in the boundary layer are largely affected by 
the thickness of the boundary layer, so the difference in the thickness of the boundary layer from the leeward side to 
the windward side should be noted. Since the disturbance frequency is inversely proportional to the boundary layer 
thickness, ^ ^ larger disturbance frequencies would be expected on the windward side of the vehicle. 


w oc 


1 

6 


(4) 


0.250 - 
0.200 - 
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This methodology, however, has limitations. As previously mentioned, the flow field around a cone at angle of 
attack is very three-dimensional. Approximating the three-dimensional flow field with a two-dimensional centerline 
slice ignores many effects including curvature, crossflow, and other three-dimensional effects that have significant 
effects on the transition behavior of a boundary layer. As a result, the analysis presented in this paper assumes that 
a two-dimensional slice along the plane of symmetry of a three-dimensional solution provides a representative flow 
field for stability analysis by minimizing the three-dimensional effects. This assumption will introduce errors into 
the analysis; the quantification of which may be addressed in future work. 

Transition onset prediction using the N factor method described is semi-empirical. The maximum N factor curve 
is computed for each condition and the N factor at transition must be determined by comparison to the measured 
transition onset location from experiment. The transition onset location for each of the ground test conditions from 
the CUBRC test series was computed using the method of zero intermittency. This method uses the time averaged 
heat transfer data from the discrete thin film gauges located on the test article to extrapolate a trend line through 
the transitional data to the point of zero intermittency, or the point where this line intersects with the laminar CFD 
prediction. This method was chosen because time traces of the thin film gauges for the CUBRC ground tests, 
which can provide more isight into the onset of transition, were not available. The uncertainty associated with the 
streamwise location of transition onset was the distance between the last gauge showing laminar behavior and the 
first gauge to deviate from the laminar prediction and exhibit transitional behavior. 

Analyzing the transition onset location from the LaRC relied more on engineering judgement because the heat 
transfer data obtained from this test was spatially continuous. The point of transition onset was chosen where the 
heat transfer data began to change slope and increase above the laminar prediction. 

Figure 5 shows how each of these methods was used to select the best estimate transition onset location and apply 
a range of transition onset uncertainty. The transition onset range shown in these figures gives the uncertainty in the 
streamwise direction along the body associated with the transition onset location described previously. It should 
be noted that the determination of the transition onset location and the associated uncertainties is very subjective. 
This publication presents only the author’s analysis. An independent analysis would likely produce different N 
factors based on transition location selection, which may or may not have an impact on the transition onset location 
predicted by the method. 




(a) CUBRC Discrete Measurements with Zero Intermittency Line, (b) LaRC Continuous Measurements with Observed Transition On- 
set. 

Figure 5. Transition onset determination methods. 
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III. Results 


Conditions were analyzed for the HIFiRE-1 geometry at ground test conditions from the CUBRC LENS I test 
facility as well as the Langley Research Center (LaRC) 20-inch Mach 6 Tunnel. The conditions analyzed can 
be found in Table 1, where the two run numbers from the LaRC tests denote the leeward and windward runs, 
respectively, since the temperature sensitive paint data at that facility was gathered only on one side of the vehicle 
for each run of this test. Lor all LaRC test data analyzed in this publication, the 1.19 mm nose radius was used. 

Lor each run, the transition onset location obtained from the experimental data was compared to the maximum 
N factor curve computed by the boundary layer stability analysis on the centerline slice of the solution. The N factor 
was extracted at three streamwise locations for each run on the windward and leeward sides: the best estimated 
transition onset location and the associated positive and negative uncertainties. Reg/Me was also extracted at the 
same locations for comparison of the empirical transition prediction methodology to the physics-based N factor 
method. Trends of the transition onset location values of N factor and Reg /Me with freestream Reynolds number 
and angle of attack were examined. 

A. Leeward Analysis 

Analysis of the boundary layer on the leeward side of the vehicle showed that it was inherently less stable than the 
windside boundary layer. Llow structures developed in the mean flow solutions that greatly affected the boundary 
layer stability. The flow did not separate, but the angle of attack effects included a ’mushroom’ of fluid build up 
along the centerline in the spanwise direction that caused inflection points within the velocity profile in the boundary 
layer. Low momentum fluid from low within the boundary layer was brought up into the higher speed flow outside 
the boundary layer and higher momentum fluid was pushed down towards the surface. The ’mushroom’ and the 
resulting inflection points in the boundary layer can be seen in the 3D mean flow solution in Ligure 6. Inflection 
points in the boundary layer are a strong instability mechanism that cause rapid transition to turbulence. Data showed 
that transition occurred shortly downstream of the location where inflection points in the boundary layer were first 
observed. 




(a) Low momentum ’mushroom’ developing in the spanwise direc- (b) Inflection points in the boundary layer profile along the centerline 
tion on the lee side of the vehicle. in the streamwise direction. 


Figure 6. Leeward Instability Mechanisms 


The location of transition onset observed during ground testing was compared to the Reg /Me computed in the 
mean flow simulation as well as the maximum N factor curve generated by STABL for each condition. The values 
of Reg /Me compared with the test data for CUBRC Run 17 and LaRC Run 19 on the leeward side of the vehicle 
can be seen in Ligure 7. Similarly, the N factor comparison for the same runs can be seen in Ligure 8. In both of the 
comparisons for the LaRC test data, there is a discrepancy between the heat transfer data obtained using the thermal 
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phosphor paint and the laminar heat transfer computed in the mean flow solution. This is due to the data resolution of 
the thermal phosphor paint which was not able to capture the increased heating on the centerline due to the vortices 
rolling up from the windward side of the vehicle. In these cases, the transition onset location was selected where the 
heat transfer data began to deviate from the expected heat transfer behavior without regard to the vortices, similar to 
the laminar heat transfer profile computed by DPLR shown in Figure 7a. 




(a) CUBRC Run 17 Leeward Data (b) LaRC Run 19 Leeward Data 

Figure 7. Ree /Me at Transition Onset 




(a) CUBRC Run 17 Leeward Data (b) LaRC Run 19 Leeward Data 

Figure 8. N Factor at Transition Onset 


The value of Reg /Mg and N factor at the transition onset location determined from the data on the leeward 
side was complied for all ground test runs to examine trends and assess uncertainties associated with each method. 
Trends with freestream Reynolds number and angle of attack were examined, but results showed only trends with 
angle of attack on the leeward side. Two runs from the CUBRC data set at 5° angle of attack were excluded from 
this analysis because the gauge resolution on the leeward side of the vehicle was not high enough to capture any 
laminar data before transition onset. This made determining the transition onset location impossible based on the 
data available using the methods described here. The angle of attack trends can be seen in Figure 9. 
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(a) 7?ee /Me Method (b) N Factor Method 

Figure 9. Correlations of leeward data obtained at transition onset 

Figure 9a shows the computed Reg /Mg at transition onset as a function of angle of attack. This empirical 
correlation methodology shows a strong trend with angle of attack, but wide scatter in the data. The mean value of 
the data obtained at transition onset results in an value of only 0.51, which emphasizes this scatter. The scatter in 
the data also increases with the angle of attack so uncertainty also increases with angle of attack. Ignoring this fact 
and placing l-a uncertainties on the data results in a value that is 82% of the mean. These l-cr uncertainties do not 
cover the spread of the data and bounding uncertainties result in an asymmetric uncertainty distribution. This may 
mean that the uncertainty in the data does not follow a Gaussian distribution, but further analysis of the uncertainties 
is needed. 

A similar correlation developed using the N factor at transition significantly reduces the uncertainties associated 
with predicting the transition using this method. Figure 9b also shows a strong trend of the N factor at transition 
onset with angle of attack but the scatter in the data is much less than that seen in the Reg /Me comparison. This 
can be seen by a larger R? of 0.83 on the mean of the data and the l-cr uncertainties on the mean span most of the 
data obtained at transition onset as well as some of the uncertainties associated with determining the transition onset 
location. The l-cr uncertainty is only 9% of the mean value for the N factor method, and a bounding uncertainty 
placed on all data is more symmetric than the Reg/M^ correlation which could mean the uncertainties in this case 
follow a Gaussian distribution more closely and be better predicted. 

B. Windward Analysis 

Ground test data on the windward side of the vehicle was compared to the maximum N factor curve and the value of 
Reg /Me and these values were extracted at the observed transition onset location and at the endpoints of the positive 
and negative uncertainties as in the leeward data analysis. The windward data showed higher heat transfer rates than 
the leeward side, but no obvious three-dimensional effects such as the low-momentum ’mushroom’ or boundary 
layer inflection points. The windward ray of gauges in the CURBC test had less gauge resolution than the ray of 
gauges on the leeward side of the vehicle which causes greater uncertainty in the transition onset location on the 
windward side. In many cases in the LaRC tests, the flow remained laminar for the entire length of the body. This 
phenomenon is counterintuitive, but has been seen in many investigations of transition on cones at angle of attack 
and is attributed to three-dimensional effects in the flow. Figure 10 shows the windward data for Run 1 1 from the 
CUBRC test series and Run 19 from the LaRC test. The maximum N factor curves for the LaRC data were also 
much lower in magnitude than the maximum N factor curves computed for the CUBRC data. 
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(a) CUBRC Run 1 1 


(b) LaRC Run 19 


Figure 10. Windward Data Comparisons 

No trends of N factor or Ree/Me with freestream Reynolds number or angle of attack were observed, so the 
uncertainties assessed on the windward side of the vehicle were examined without regard for a physical flow or 
geometric parameter. The mean value of the data extracted at the transition onset location and 1-cj and bounding un- 
certainties were examined. Comparisons for the N factor at transition onset for the LaRC data and the CUBRC data 
differed by a factor of at least four based on the two LaRC runs which exhibited transitional behavior on the wind- 
ward side. As a result, the LaRC data and CUBRC data were analyzed separately. More than two transitional data 
points from the LaRC test series should be examined to create any correlations and draw meaningful conclusions. 
As a result, only data from the CUBRC test is presented here. 

The results of the Ree/Me empirical correlation for the CUBRC data, seen in Figure 11a, shows more scatter 
in the data and larger uncertainties than those observed in the evaluation of the N factors computed at transition 
(Figure 1 lb). The l-u uncertainty in the Rcq /M^ methodology is reduced when using the N factor method. Opposite 
of the findings on the leeward side of the vehicle, the bounding uncertainties in the Ree /Me empirical methodology 
are more symmetric than those for the N factor methodology. A more detailed analysis of the distribution of the 
uncertainties around the mean value will be required to address this finding. 
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(a) Reel Me Method 


(b) N Factor Method 


Figure 11. Correlations of windward data obtained at transition onset from CUBRC data. 
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For the LaRC windward data, two runs showed transitional behavior - Run 30 at 5° angle of attaek and Run 65 
at 3° angle of attaek. The N faetors at transition were 1.5 and 0.6, respeetively, whieh were signifieantly less than 
the N faetors at transition eomputed from the CUBRC eonditions. Additionally, for all of the other runs in the LaRC 
data set that did not show transition in the heat transfer data on the windward side, the maximum N faetor eurves did 
not reaeh above 2. The differenee in the relative nose radius between the CUBRC and LaRC test artieles previously 
diseussed eould explain this differenee in part, but many other faetors eould be eontributing to the large differenee 
in N faetors. 


IV. Conclusion 

Ground test data for the HIFiRE-1 geometry from the CUBRC LENS I refleeted shoek tunnel and Eangley 
Researeh Center 20-ineh Maeh 6 Tunnel was used to examine the boundary layer stability and transition onset 
eharaeteristies at non-zero angles of attaek. Diserete heat transfer measurements from thin film gauges from the 
CUBRC data sets and spatially eontinuous heat transfer measurements from thermal phosphor paint during the 
EaRC tests were used to determine transition onset loeation. The transition onset loeation was used to extraet 
the boundary layer property Reg/Mf,, whieh is eommonly used to prediet boundary layer transition in empirieal 
eorrelations, from the laminar mean flow solution and the maximum N faetor eomputed from the boundary layer 
stability analysis. Windward and leeward sides of the vehiele were analyzed separately. 

On the leeward side of the vehiele, a strong trend of both transition onset indieation parameters with angle of 
attaek was observed. This eorrelates strongly to the development of infleetion points in the boundary layer that 
are a strong instability meehanism and lead to rapid breakdown to turbulenee. The windward side of the vehiele 
transitioned farther down the body than on the leeward side, but no trends were observed with angle of attaek or 
freestream Reynolds number. Using both methods of predieting boundary layer transition onset, uneertainty analysis 
was performed. The N faetor method in all eases presented here showed a stronger eorrelation with less seatter in the 
data and smaller uneertainties than the Reg/Mf, method. However, bounding uneertainties were asymmetrie about 
the mean in some eases, and a more detailed analysis that addresses the distribution of the uneertainty around the 
data must be performed. 

In general, the maximum N faetor eurves eomputed for the EaRC test data had lower values of N faetor than 
the eurves generated for the CUBRC test data. On the windward side of the vehiele, the N faetors between the 
two faeilities differed by at least a faetor of four. This eould be due to the differenee in the relative nose radius of 
the vehieles. The EaRC test artiele had a larger nose radius relative to the length of the vehiele than the CUBRC 
test artiele and the flight vehiele. A larger nose radius has been shown to stabilize the boundary layer and delay 
transition, so this eould be eontributing to the differenees in the eomputed N faetors. Eurther analysis would be 
required to determine the eause of the diserepaney with more eertainty. 

The results of this investigation showed a deerease in uneertainty when using a physies-based transition pre- 
dietion methodology when eompared to an empirieal eorrelation based on Reg/Me for the HIEiRE-1 slender eone 
geometry at ground test eonditions. The analysis at non-zero angles of attaek presented here eould be expanded in 
the future to ineorporate the results of the HIEiRE-1 flight test as most of the transition data obtained during the 
flight was at non-zero angles of attaek. Ineluding the flight data in future analysis would also begin to address the 
issue of flight traeeability in boundary layer transition onset predietions. 
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